Electromagnetic waves destabilized by runaway electrons in near-critical electric fields

Runaway electron distributions are strongly anisotropic in velocity space. This anisotropy is a source of free energy that may destabilize electromagnetic waves through a resonant interaction between the waves and the energetic electrons. In this work we investigate the high-frequency electromagnetic waves that are destabilized by runaway electron beams when the electric field is close to the critical field for runaway acceleration. Using a runaway electron distribution appropriate for the near-critical case we calculate the linear instability growth rate of these waves and conclude that the obliquely propagating whistler waves are most unstable. We show that the frequencies, wave numbers and propagation angles of the most unstable waves depend strongly on the magnetic field. Taking into account collisional and convective damping of the waves, we determine the number density of runaways that is required to destabilize the waves and show its parametric dependences.Comment: 22 pages, 11 figures, to be published in Physics of Plasma

Reducing systematic errors in time-frequency resolved mode number analysis

The present paper describes the effect of magnetic pick-up coil transfer functions on mode number analysis in magnetically confined fusion plasmas. Magnetic probes mounted inside the vacuum chamber are widely used to characterize the mode structure of magnetohydrodynamic modes, as, due to their relative simplicity and compact nature, several coils can be distributed over the vessel. Phase differences between the transfer functions of different magnetic pick-up coils lead to systematic errors in time- and frequency resolved mode number analysis. This paper presents the first in-situ, end-to-end calibration of a magnetic pick-up coil system which was carried out by using an in-vessel driving coil on ASDEX Upgrade. The effect of the phase differences in the pick-up coil transfer functions is most significant in the 50-250 kHz frequency range, where the relative phase shift between the different probes can be up to 1 radian (~60{\deg}). By applying a correction based on the transfer functions we found smaller residuals of mode number fitting in the considered discharges. In most cases an order of magnitude improvement was observed in the residuals of the mode number fits, which could open the way to investigate weaker electromagnetic oscillations with even high mode numbers

Quasi-linear analysis of the extraordinary electron wave destabilized by runaway electrons

Runaway electrons with strongly anisotropic distributions present in post-disruption tokamak plasmas can destabilize the extraordinary electron (EXEL) wave. The present work investigates the dynamics of the quasi-linear evolution of the EXEL instability for a range of different plasma parameters using a model runaway distribution function valid for highly relativistic runaway electron beams produced primarily by the avalanche process. Simulations show a rapid pitch-angle scattering of the runaway electrons in the high energy tail on the $100-1000\;\rm \mu s$ time scale. Due to the wave-particle interaction, a modification to the synchrotron radiation spectrum emitted by the runaway electron population is foreseen, exposing a possible experimental detection method for such an interaction

Experimental investigation of the radial structure of energetic particle driven modes

Alfv\'en eigenmodes (AEs) and energetic particle modes (EPMs) are often excited by energetic particles (EPs) in tokamak plasmas. One of the main open questions concerning EP driven instabilities is the non-linear evolution of the mode structure. The aim of the present paper is to investigate the properties of beta-induced AEs (BAEs) and EP driven geodesic acoustic modes (EGAMs) observed in the ramp-up phase of off-axis NBI heated ASDEX Upgrade (AUG) discharges. This paper focuses on the changes in the mode structure of BAEs/EGAMs during the non-linear chirping phase. Our investigation has shown that in case of the observed down-chirping BAEs the changes in the radial structure are smaller than the uncertainty of our measurement. This behaviour is most probably the consequence of that BAEs are normal modes, thus their radial structure strongly depends on the background plasma parameters rather than on the EP distribution. In the case of rapidly upward chirping EGAMs the analysis consistently shows shrinkage of the mode structure. The proposed explanation is that the resonance in the velocity space moves towards more passing particles which have narrower orbit widths.Comment: submitted to Nuclear Fusio

Special behavior of alkali beam emission spectroscopy in low-ion-temperature plasma

Beam emission spectroscopy (BES) is a powerful plasma diagnostic method especially suited for the measurement of plasma density and its fluctuations. As such, synthetic BES codes are regularly used to aid the design or utilization of these diagnostic systems. However, synthetic diagnostics can also be used to study the method in previously not yet explored operational conditions. This paper presents such an analysis utilizing the RENATE-OD synthetic diagnostic code for a hypothetical alkali BES system on the HSX stellarator. HSX is a device featuring an unusual operating regime in the world of fusion devices due to the low ion temperature and low plasma density. It was found that BES shows unusual tendencies in these conditions. The relation between beam energy and plasma penetration in low-ion-temperature plasma, together with unique emission features facilitated by low-density plasma, and the underlying reasons behind these features are explored in this paper

Validity of models for Dreicer generation of runaway electrons in dynamic scenarios

Runaway electron modelling efforts are motivated by the risk these energetic particles pose to large fusion devices. The sophisticated kinetic models can capture most features of the runaway electron generation but have high computational costs which can be avoided by using computationally cheaper reduced kinetic codes. In this paper, we compare the reduced kinetic and kinetic models to determine when the former solvers, based on analytical calculations assuming quasi-stationarity, can be used. The Dreicer generation rate is calculated by two different solvers in parallel in a workflow developed in the European Integrated Modelling framework, and this is complemented by calculations of a third code that is not yet integrated into the framework. Runaway Fluid, a reduced kinetic code, NORSE, a kinetic code using non-linear collision operator, and DREAM, a linearized Fokker-Planck solver, are used to investigate the effect of a dynamic change in the electric field for different plasma scenarios spanning across the whole tokamak-relevant range. We find that on time scales shorter than or comparable to the electron collision time at the critical velocity for runaway electron generation kinetic effects not captured by reduced kinetic models play an important role. This characteristic time scale is easy to calculate and can reliably be used to determine whether there is a need for kinetic modelling, or cheaper reduced kinetic codes are expected to deliver sufficiently accurate results. This criterion can be automated, and thus it can be of great benefit for the comprehensive self-consistent modelling frameworks that are attempting to simulate complex events such as tokamak start-up or disruptions

Spatiotemporal analysis of the runaway distribution function from synchrotron images in an ASDEX Upgrade disruption

Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analysed using the synchrotron diagnostic tool Soft and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quench, together with an avalanche runaway tail which has an exponentially decreasing energy spectrum. We find that, although the avalanche component carries the vast majority of the current, it is the high-energy seed remnant that dominates synchrotron emission. With insights from the fluid-kinetic simulations, an analytic model for the evolution of the runaway seed component is developed and used to reconstruct the radial density profile of the RE beam. The analysis shows that the observed change of the synchrotron pattern from circular to crescent shape is caused by a rapid redistribution of the radial profile of the runaway density